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Abstract 

We investigate the 0{as) corrections to e+e^ — ^ J/tpl^il)') + Xcj {J = 0,1,2) in the NRQCD 
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pq ' factorization approach. These next-to-leading order (NLO) corrections are calculated at the level 

of helicity amplitude. We have made a detailed analysis for both polarized and unpolarized cross 

,J^ ■ sections, and compared our predictions with the measurements at the B factories. We also de- 

^~^, rive the asymptotic expressions for each of the NLO helicity amplitudes, and confirm the earlier 

0-i| speculation that at NLO in ct^, the double logarithm of type In {s/rn^) appearing in the NRQCD 

(— I \ short-distance coefficient is always associated with the helicity-suppressed channels. 
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I. INTRODUCTION 

The studies of hard exclusive reactions have historically played an important role in the 
development of perturbative Quantum Chromodynamics (pQCD). The standard theoretical 
tool is the light-cone (coUinear factorization) approach, which is based on the expansion 
in powers of 1/Q, where Q^ characterizes the hard momentum transfer [l|, |2|. The classic 
applications of light-cone factorization are exemplified by the 71—7 transition form factor, tt 
electromagnetic (EM) form factor [l|, |2[, and the B meson exclusive decays |3|, ^ , to which 
a vast amount of literature has been devoted. 

Recent advancement of the high-luminosity high-energy collider facilities renders the 
study of hard exclusive processes involving heavy quarkonium a fertile new research frontier. 
Perhaps a great amount of interest toward this topic was triggered by the observation of 
several double-charmonium production processes in two B factories several years ago p, |6|. 

In addition to the light-cone approach, there also exists another relatively new factor- 
ization approach, the NRQCD factorization [7|, which is tailor-made to tackle quarkonium 
production and decay processes. Explicitly exploiting the nonrelativistic nature of quarko- 
nium, the NRQCD factorization approach allows one to express the amplitude of an exclusive 
quarkonium production reaction in terms of an infinite sum of products of short- distance 
coefficients and the vacuum-to-quarkonium NRQCD matrix elements, whose importance is 
organized by the powers of f , the typical velocity of a heavy quark inside the quarkonium. 

Perhaps the most famous example of double charmonium production is e^e~ — )■ J/tp + r)c, 
with the interest originally spurred by the alarming discrepancy between the Belle mea- 



surement [5| and the leading-order (LO) NRQCD predictions |8HlO|. In the following years, 
this process has been intensively studied in both of the NRQCD and light-cone frame- 
works [sHis . 



One of the important steps toward alleviating the discrepancy between data and theory for 
e~^e~ — )■ J/ip+r]c is the discovery of significant and positive NLO perturbative corrections [llj. 



12| . This NLO calculation was performed in NRQCD factorization. In contrast, due to some 
long-standing theoretical difficulty inherent to this helicity-suppressed process, by far no one 
has successfully worked out the NLO correction to this process in the light-cone approach. 
Besides e+e^ — )■ J/tp + rjc, two B factories have also measured several other double- 
charmonium production processes p, |6|, notably e"*"e~ — )■ J/ip + Xcj (J = 0, 1, 2). Recently 



in Ref. [16|, the NLO perturbative correction was performed for e~^e~ — )■ J/ip + Xco at LO in 
V and a rather large K factor was reported to bring the theory prediction in better agreement 
with the measurement. Nevertheless, at this stage, both B factory experiments have not 
clearly observed any J/ip + Xci,2 events yet, even with latest data set [17], and only an upper 
bound for the production cross sections is placed. 

The purpose of this work is to carry out a comprehensive next-to-leading order (NLO) 
perturbative analysis to e"^e~ — > J/ip^ip') + Xcj {J = 0, 1, 2) at the lowest order in v. We 
will work out the NLO corrections to each of the helicity amplitudes, so that we can address 
their impact on both the polarized and unpolarized cross sections. 

We find a large positive NLO perturbative correction to e^e~ — >■ J /ip + XcOi which is 
helpful to bring the predicted cross section closer to the B factory measurements. However, 
our NLO predictions to ip' + Xco cross section seems still significantly below the central value 
of the Belle measurement. 

On the other hand, the impact of NLO corrections to e^e~ — )■ J /ip + Xci,2 seems to 
be rather modest, even with their signs uncertain. Our predicted cross sections for these 



processes are about 5 or 6 times smaller than that for e^e~ -4 J/ip + Xco- Hopefully the 
future Super B factory, with much higher luminosity, will eventually observe these two 
channels. 

Our studies of polarized cross sections reveal that the bulk of the total cross section comes 
from the (0, ±1) helicity states for e+e^ — t- J/tp + Xd, and from the (0, 0) and (±1, 0) helicity 
states for e~^e~ — )■ J/ip + Xc2- It will be interesting for the future Super B experiments to 
test these polarization patterns. 

From the theoretical perspective, we have also presented the analytic asymptotic expres- 
sions of all the ten NLO helicity amplitudes for the e'^e'^ — > J/ip + Xcj {J = 0, 1, 2) processes. 
The logarithmical scaling-violation pattern is recognized in these NLO asymptotic expres- 
sions, which supports for the speculation made in Ref. |l8|: The hard exclusive processes 
involving double-quarkonium at leading twist can only host the single collinear logarithm 
Ins /ml at one-loop, while those beginning with the higher twist contributions are always 
plagued with double logarithms of form In^ s/m^. It would be theoretically interesting to 
reproduce these asymptotic expressions by using the light-cone approach. 

The rest of the paper is organized as follows. In Section [TTl we present a concise review 
on the helicity amplitude formalism and helicity selection rule, taking the process e~e~^ — )■ 
J/ip + Xcj {J = 0, 1, 2) as the example. In Section UTTt we rederive the tree-level predictions 
in the NRQCD factorization approach, and present the LO expressions for all the involved 
helicity amplitudes. In Section llVt we first elaborate on some key technical issues about 
the NLO perturbative calculations, then present the asymptotic expressions for the NLO 
corrections for all the encountered helicity amplitudes. The pattern of the logarithmic scaling 
violation is recognized, and their implication with the light-cone approach is discussed. We 
devote Section |V] to explore the phenomenological impact of our NLO predictions on the B 
factory measurements. A comprehensive numerical study is performed for both unpolarized 
and polarized cross sections of the e'e"*" — )■ J/ip + Xcj {J = 0, 1,2) processes. Finally, we 
summarize in Section |Vl 

II. POLARIZED CROSS SECTIONS AND HELICITY SELECTION RULE 

It is often desirable to know the polarization information for an exclusive reaction, es- 
pecially for the double-charmonium production process considered in this work. To fulfill 
this goal, one can employ the helicity amplitude formalism [l9|, which appears to be a quite 
convenient and advantageous tool both experimentally and theoretically. 

Suppose the e~ and e"*" beams are aligned along the z direction, bearing the invariant mass 
of y^. We will choose to work in their center-of-mass frame. Let 6 denote the angle between 
the moving directions of J/ip and e~, and |P| signify the magnitude of the momentum 
carried by the J/ip (Xcj) [Note that |P| = As (s, Mj, ,, M^^ J/(2^), where X{x,y,z) = 
x^ +y'^ + z^ — 2xy — 2yz — 2zx] . Assume the outgoing J/ip and Xcj carry definite helicities 
of Ai, A2, respectively. First let us imagine the process of a virtual photon decay into J/ip 



and Xcj- The differential decay rate can be expressed as [19|, |20 



dr[r(g.)^J/^(Ai)+Xcj(A2)] _ |P| I 1 ..,2, , ,2 

where A = Ai — A2, ^2 is the spin projection of the virtual photon along the z axis. A'l^ X2 
is the desired helicity amplitude, which encodes all the nontrivial dynamics. The angular 



distribution is fully dictated by the quantum numbers Sz and A through the Wigner rotation 
matrix dl^^,{6). Note angular momentum conservation constrains that |A| < 1. 

The number of independent helicity amplitudes can be greatly reduced, thanks to the 
parity invariance: 

Km = (-l)'^-A„-A.- (2) 

Consequently, the two helicity amplitudes related by simultaneously flipping the helicities 
of J/tp and Xcj bear the equal magnitude. 

It is straightforward to convert the differential decay rate ([1]) to the corresponding pro- 
duction cross section in e^e~ annihilation: 

da[e+e- ^ J/tpjXi) + Xcj(A2)] _ 27ra y^ dT[YiSz) ^ J /^jX,) + Xcj iX2)] 

5z=±l 

"SsHv^y"^^-^^' ""[ sin^^, (A = 0) ^^^ 

where the polarizations of e~ and e"*" have been averaged over. Note that we only need to 
sum over two transverse polarizations of the virtual photon, as guaranteed by the helicity 
conservation in QED for a massless electron. This selective summation is the cause for the 
anisotropic angular distribution patterns in ([3]). 

It is now straightforward to acquire the unpolarized cross sections. Integrating ([3]) over 
the polar angle 6 and including all the allowed helicity states, one finally arrives at: 

^[J/^ + Xco] = £^ (^) {\aW + 2 \A%\') , (4a) 

cr[J/^ + Xci] = £-, (^) (2 1^1,0 r + 2 1^0,1 r + 2 I^Ur) , (4b) 

cr[J/^ + Xc2] = £-, (^) {\Al/ + 2 \Al/ + 2 \Al,\' + 2 \Al\' + 2 \Al/) .(4c) 

There are 2, 3 and 5 independent helicity amplitudes for 7* — )■ J/ip + Xcj {J = 0, 1,2), 
respectively, as required by the angular momentum conservation. Note that parity trans- 
formation property (|2]) forbids the appearance of the AIq amplitude. In equation @, we 
include a factor of 2 to account for the parity-doublet contribution. 

Aside from its great technical usefulness, the helicity amplitude formalism can also shed 
important light on the dynamics underlying exclusive reactions. In particular, each helicity 
amplitude possesses a definite power-law scahng in 1/s, controlled by the helicity selection 
rule [21]. At asymptotically large ^/s, the polarized cross section for our double- charmonium 
production process scales as [8|: 



m 



2\ 2+IA1+A2I 



a[JMAi) + Xcj(A2)] ~ aV(^J , (5) 

here v denotes the characteristic velocity of charm quark inside a charmonium. 

Equation ([5]), combined with angular momentum conservation, implies that the helicity 
state which exhibits the slowest asymptotic decrease, i.e. a oc 1/s^, is the one that conserves 
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FIG. 1: One sample LO diagram and five sample NLO diagrams that contribute to 7* — t- J/ijj + XcJ- 



the hadron helicities |Ai + A2I = 0, or equivalently, (Ai, A2) = (0,0). For each unit of the 
violation of this law, there is a further suppression factor of 1/s. At sufficiently high energy, 
only the (0, 0) helicity state perhaps needs to be retained for phenomenological purposes. 
Note that in NRQCD factorization language, the charm quark is treated as heavy, and in 
fact its mass acts as the agent of violating the hadron helicity conservation. 

Once beyond the lowest order in a^, the power-law scaling specified in ([5]) will in general 
receive mild modifications due to the In s terms. This logarithmic scaling violation will be 
examined in detail in Section IIVI 



III. LO PREDICTIONS FOR POLARIZED CROSS SECTIONS 



At LO in V expansion in the NRQCD approach, one can factorize the amplitude of 
7* — >■ J /if] + Xcj as the product of the short-distance coefficients and the nonperturbative 
NRQCD matrix elements. To a good extent, these nonperturbative matrix elements are well 
approximated by the phenomenological (derivative of) wave function at the origin, i?j/^(0) 
and i?^^^(0). At LO in a^, the short- distance coefficients have been obtained [8|, Isj by 



/3q(1) 



3p(l)^ 



calculating the quark amplitude 7* — )■ cc( S\ ) + cc{ Pj ) using the covariant projection 
technique |22|, |23| . 

There are in total 4 Feynman diagrams at LO in as, one of which is depicted in Fig. [TJi). 
For simplicity, we have neglected the QED fragmentation diagrams, while their effect is 
modest. To expedite the extraction of the corresponding helicity amplitudes, we have con- 
structed 10 helicity projection operators. It is convenient to parameterize the LO helicity 
amplitude as 



.jio) _ ^ee,asCFRj/4^)R'^JO) Ui+\x^+X2\) ,J (^^ 

•^ Ai,A2 - ZI^ ^ (^XiM^h 



(6) 



where Cc = I is the electric charge of the charm quark, and for brevity, we have introduced 



a dimensionless variable r: 

Am 



r ^ —^. (7) 



To make the scaling behavior in ([5]) transparent, we have explicitly stripped off a factor 



f2 



i(l+|Ai+A2|) 



in ([6]), so that the reduced hehcity amplitude c^^j^^l'^); ^ dimensionless function, 
will start with a 0{1) constant. Concretely, these c'l_^ a2('") functions are 

0° (,(r) = -1 - lOr + 12r2 c%{r) = -9 + 14r, (8a) 

ci_^(r) = 76(2 - 7r) cj o(r) = -V6r cl-^{r) = 2^6(1 - 3r), (8b) 

4o(r) = y2(-l + 2r + 12r2) ^^^^(r) = V6(-l + 5r) ^^(r) = v^(ll - 3r), 
c^_i(r) = 2^6(1 - 3r) c^ 2(r) = -2^3. (8c) 



These 10 helicity amplitudes agree with those given in Ref. [8|. It is interesting to note 
that the tree-level J /tp{l) + Xci(O) amplitude constitutes an exception in that it accidently 
receives an extra suppression factor than expected from the helicity selection rule. 

Plugging ([6]) into (jl]), we obtain the LO NRQCD predictions for the polarized cross 
sections: 

a(°)[j/^(Ao +xc.(A.)] = ''i';!?"' (^) R%,miM i<,A.wr ■ (9) 

IV. NLO PERTURBATIVE CORRECTIONS TO POLARIZED CROSS SEC- 
TIONS 

In this section, we first sketch some technical issues about the NLO perturbative calcula- 
tions, then present the asymptotic expressions of the NLO corrections for all the encountered 
helicity amplitudes. 



A. Description of the calculation 



We first employ the Mathematica package FeynArts [2J, |25| to generate Feynman 
diagrams and amplitudes for the parton process 7* — )■ cc(Pi)+cc(P2) to NLO in a^. Feynman 
gauge is used throughout this calculation. In total there are 20 two-point, 20 three-point, 
18 four-point, and 6 five-point one-loop diagrams, some of which have been illustrated in 
Fig. [H We then apply the covariant projector technique [22, l23| to enforce that two cc pairs 



form the spin-triplet, color-singlet states, with the Dirac and color traces handled by the 
package FeynCalc [26]. In the next step, we proceed to expand the integrand in powers of 
quark relative momenta, g, to project out the leading S-wave and P-wave orbital- angular- 
momentum contributions. 

Making expansion in q prior to carrying out the loop integration, in the spirit of method 
of region 27|, amounts to directly deducing the NRQCD short-distance coefficients, i.e., the 



contributions solely stemming from the hard region (/c^ > rn^). In practice, this procedure 
is far simpler than the conventional perturbative matching procedure, since by this way we 
will no longer be distracted by the effects from the low-energy regions, e.g., from the soft 



(A;^ ~ mv), or potential {k^ ~ mv'^, |k| ~ mv) regions. Consequently, with this strategy, one 
will never confront Coulomb singularity in the one-loop integrals. 

In Ref. J28[], an all-order- in-a^ proof for exclusive quarkonium production has been out- 
lined in the NRQCD factorization context. It is argued that, at lowest order in v and to all 
orders in a^, NRQCD factorization holds for S'-wave plus P-wave quarkonia production in 
e"'"e~ annihilation. In particular, this implies that the NRQCD short-distance coefficients 
associated with e^e~ — )■ J/ip + Xcj should be free from any infrared singularity at NLO in 
as ^. Indeed, an earher NLO perturbative calculation for e'^e" — )■ J/ip + Xco has explicitly 
verified this pattern [l6|. 

At this stage, we apply the specifically-designed helicity projection operators to project 
out 10 corresponding helicity amplitudes. This operation brings forth great simplification, 
for all the polarization vectors (tensors) have been eliminated, and the numerators in loop 
integrals now become the Lorentz scalars comprised entirely of external and loop momenta. 
Note these integrals in general contain propagators with quadratic power due to the projec- 
tion of the P-wave state. 

To proceed, we use the Mathematica package FIRE |32| and our self-written Math- 
EMATICA code to reduce the general higher-point one-loop integrals into a set of Master 
Integrals (MI). Fortunately, as a great bonus of having Taylor-expanded the integrand in 
powers of q, together with having applied the helicity projectors, there are only three linearly- 
independent Lorentz scalars in loop integrals: /^, / ■ Pi and / ■ P2, where /, Pi, P2 stand for 
the loop momentum, the momenta carried by the cc(^S'| ) pair and by the cc{^Pj ) pair, 
respectively. Thanks to the integration-by-part algorithm built in FIRE, it turns out that 
all the involved MI become just the usual 2-point and 3-point scalar integrals. All the en- 



countered 3-point scalar integrals can be found in Appendix of Ref. [12| . We have used the 
package LoopTools [SS^ to numerically check the correctness of these integrals. 

Throughout this calculation we adopt dimensional regularization to regularize both the 
ultraviolet (UV) and infrared (IR) singularities. In accordance with the LSZ reduction 
formula, one needs to multiply the tree-level amplitude by (a/Z^)'^, where Z^ denotes the 
residue of the charm quark propagator near its pole. This contribution is represented by 
Fig. [1]/). In Feynman gauge, the residue is given by 



Zc — I — Cp-r- 

4lT 



2 



2 , inn 

+ 37ij + 3 In — ^ + 4 



euv eiR rn'^c 



+ Oiai), (10) 



where je is the Euler constant, and Cp 



2Nc ■ 

In addition, we also need to replace the bare charm quark mass and the bare QCD 



^ An uncanceled infrared singularity has been discovered in NLO calculation for the B — ^ Kxcj 2^ |30| , 
which has been interpreted as a failure of NRQCD factorization. This symptom hints that the lower-energy 
degrees of freedom, the ultrasoft (k ~ mv'^) region, has been erroneously missed, and one should invoke 
the even lower energy effective theory, the potential NRQCD, to remedy this problem. Once the ultrasoft 
mode from higher Fock state |cc(^5'} )g) is consistently taken into account, this infrared divergence can 
be tamed 3l|. 



coupling constant in the tree-level amplitude by 
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bare 
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47r V^uv rnt. 3 
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MS 



A 
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- 7s + ln(47r) + 0{a 



(11a) 
(lib) 



where (3q = yA^c — 3^/ is the one-loop coefficient of the QCD (3 function, and Uf = 4 denotes 
the number of active quark flavors. 

When adding the contributions from all the diagrams, we find that the ultimate NLO 
expressions for each of the 10 helicity amplitudes are UV and IR finite. 



B. Analytic expressions of NLO helicity amplitudes 

For clarity, we parameterize the NLO helicity amplitude as follows: 



A 



j(i) 

Ai,A2 



«., 



— K 



TT 



Ai,A2 



rA\A 



,j(o) 

Ai,A2' 



(12) 



where the LO helicity amplitude A^^^ is defined in 1^. The dimensionless quantity K^^ ^^, a 
function of r and the renormalization scale /i, can be regarded as the reduced NLO helicity 
amplitude. It necessarily encompasses all the loop dynamics. Since helicity selection rule 
has been tacitly embodied in A^^^ in f lT2|) . we expect that the K functions will exhibit slower 
asymptotic decrease than any power-law scaling in r. 

The reduced NLO helicity amplitudes K are in general complex-valued. Their analytic 
expressions are somewhat lengthy and will not be reproduced here. On the other hand, in 
Figs. 121 [31 and IH we explicitly display their profiles over a wide range of r for each of the 
10 helicity amplitudes. 

Combining ([6]), ( IT2l) and (^, we can deduce the NLO perturbative correction to the 
polarized cross section: 



a«[JMAi) + Xcj(A2)] = 2(^)Re| 






^ \ I ^(0) 



'[J/^(Ai) + Xcj(A2)], (13) 



where a^^^ is given in Q. To the desired NLO accuracy, the imaginary part of the K function 
does not contribute. We now have improved prediction o"nlo = (^^^^ + cr^^^- 

It is theoretically curious to know the asymptotic behaviors of these reduced helicity 
amplitudes in the limit y/s ^ 171^. As was mentioned before, we anticipate to see the 
logarithmic scaling violation at NLO to the naive power-law scaling rule (|5|). Moreover, 
performing the asymptotic expansion in NRQCD short-distance coefficients is theoretically 
very appealing, since it is equivalent to disentangling the effects occurring at the "hard" scale 
(virtuality ~ s) from the lower-energy "coUinear/soft" effects (virtuality ~ fn^), by which one 
can intimately link the NRQCD factorization approach and the light-cone approach [l8|, |34 . 
Such asymptotic expansion has been carried out for several double-quarkonium production 
processes at one-loop level [l8|, |35|, |36| and some general patterns have been recognized. 
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FIG. 2: The profiles of two reduced NLO helicity amplitudes K^ ^^ (for 7* — )• J/ip + Xco) with ^/s 

varied in a wide range. We take /x = -^ and rric = 1.5 GeV. The solid curves correspond to the 
exact results, and the dashed curves represent the asymptotic ones as given in ()14p . The vertical 
mark is placed at the B factory energy ^/s = 10.58 GeV. 
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FIG. 3: The profiles of three reduced NLO helicity amplitudes K^^ ^^ (for 7* — )■ J/ip + Xci)- The 
asymptotic curves are taken from (J15p . The parameters are the same as in Fig. [2l 



The asymptotic expressions of two reduced NLO helicity amplitudes for J/ip + Xco are 



__L(46 + vr2-401n2 + 331n22) + ^(^/3o-y + ^ln2^, 



(14a) 
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(14b) 



The asymptotic expressions of the three reduced NLO helicity amplitudes for J/ip + Xci 
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FIG. 4: The profiles of five reduced NLO helicity amplitudes K^^ ^^ (for 7* — )• J/ip + Xc2)- The 
asymptotic curves are taken from (J16p . The parameters are the same as in Fig. [2l 
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51nV + (7-21n2)lnr- 19 + 27r^ + 751n2 - 21 In^ 2 
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(15a) 



^0,1 (^' — ) = —\ — InV- (46-991n2)lnr + /3o(61n— + 121n2+ 13 

V ^ /asym 96 I 2 \ S 



— (616 + 747r^ - 1696 ln2 + 303 In^ 2) + in{25\nr + 6/5o - 46 + 991n2) }, (15b) 



6 



-"(■4).,M 
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The asymptotic expressions of the five reduced NLO helicity amplitudes for J/ip + Xc2 



10 



are 

Kj„(4)^^^^^-l(4-ln2)ln. + ^(ln^ + ^ + 21n2 

_l(64 + 7r2 + 104 In 2 + 33 In^ 2) + '-^f^, - ^ + ^^2)' ^^^^) 

^0,1 (^' — ) = t\ — lnV-(22-431n2)lnr + /3o(31n^ + 8 + 61n2 

' V ^ J asym 48 |^ 2 \ S 

— (284 + 307r2-3801n2 + 1591n2 2)+i7r(131nr + 3/3o-14 + 431n2) J>, (16b) 
6 

KlJr,^\ =-^<21n2r + ^(5 + 81n2)lnr + /3o(^ln^ + ^ + 31n2 

-^(291-87r2 + 1711n2 + 3121n22)+i7rUlnr + ^/3o + y + ^ln2) }, (16c) 

^1,1 (^' — ) = — <41nV- (46-621n2)lnr + /3o(61n— + 13 + 121n2 
■^ v. 

— (274 + 277r2 - 316 In 2 + 9 In^ 2) + i7r(8 In r + 6/3o - 46 + 62 In 2) J> , (16d) 
^1,2 (^A) " -^<21nV + ^(l + 131n2)lnr-/3o('ln^ + ^ + 21n2 

1' -.2 , ... ...,„o , oo.,„2.. . /. . ^ 26 



+ -(-77r^ + 140-1041n2 + 2371n^2)-i7rl /3o + 3-— ln2 ) J>. (16e) 

Now we can make several interesting observations from Eqs. flMl) . f lTSj) . f lT6|) . One con- 
firms that the scahng violation is indeed of the logarithmic form. For the hadron-helicity- 
conserving channels such as J/ip{0) + Xco,2(0), the leading scaling behavior of the real part 
of the K function is governed by a single logarithm of r; for all the remaining helicity- 
suppressed channels, the leading asymptotic behaviors of the respective K functions are all 
proportional to In^ r. 



This pattern lends further support for the speculation made in Ref. [18|: The hard exclu- 
sive processes involving double-quarkonium at leading twist can only accommodate the single 
collinear logarithm at one-loop, while those beginning with the higher twist contributions 
are always plagued with double logarithms. For example, the NLO NRQCD short-distance 
coefficients of the helicity-suppressed reactions e"'"e~ — )■ J/ip + Vc and rjb — )■ J/'ipJ/ip are both 
found to contain the double-logarithm term [l8|, |36!|. 

Most of the studied double-charmonium production processes are of the helicity- 
suppressed type. In this regard, the hehcity channels e+e~ — )■ J/ipiO) + Xco,2(0) constitute 
the rather rare examples that the leading-twist contribution dominates. In such a situation, 
by resorting to the leading-twist collinear factorization theorem, one can employ the light- 
cone approach to efficiently reproduce the asymptotic expressions given in (I14ap and (I16ap . 
very much like what is achieved for Be electromagnetic form factor at NLO in ag [18|. Note 
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that the single logarithm in these two channels have identical coefficient oc 4 — In 2. This 
coefficient can be readily reconstructed in light-cone approach, with the aid of the Efreniov- 



Radyushkin-Brodsky-Lepage (ERBL) evolution equation [37|, l38||. Moreover, following the 
strategy of [SJ], by employing this evolution equation, one can systematically identify and 
resum the leading collinear logarithms in these amplitudes to all orders in a^. 

In contrast, it remains to be an open challenge for light-cone approach to reproduce 
these pro cess- dependent double logarithms appearing in NLO NRQCD short- distance co- 
efficients, which seem to have resulted from the overlap between collinear and endpoint 
singularities [l8|. Note that the end-point singularity is a long-standing problem in light- 
cone framework, which has essentially hindered our capability of performing the complete 
NLO perturbative calculation beyond leading-twist using this approach. Perhaps the ffist 
major progress is to successfully reproduce the asymptotic expressions related with those 
(0, 1) and (1,0) hehcity channels, where only twist-3 effects need be considered. 

For reader's convenience, all the asymptotic results of the reduced helicity amplitudes 
are also shown in Figs. [21 [31 IH in juxtapose with the corresponding exact NLO results. As 
can be seen clearly, for those higher-twist helicity channels, the asymptotic results tend to 
converge with the exact ones decently well at relatively lower ^/s (say, at B factory energy) 
than for the leading-twist channels. This seems to be a general feature, which has also been 
observed in several other double charmonium production processes [l8|, |35|, ISg . 



To conclude this section, we mention a peculiar phenomenon affiliated with the helicity 
channel 7* — )■ J / ip (1) + Xci{0) ■ Recall that this helicity amplitude at LO has been suppressed 
by an unwanted factor of r, as can be seen in ( JSbll . Nevertheless, as one may readily tell 
from (I15ap . at NLO in as, this helicity amplitude recovers the correct power- law scaling 
as dictated by the helicity selection rule! This implies that at very high energy, the NLO 
contribution is far more significant than the LO piece for this helicity channel, despite its 
extra suppression by a^. It might be also worth noting that, unlike all the other asymptotic 
expressions, the renormalization logarithm PqIu fi^/s is absent in fll5a ). This is very similar 
to what is found in NLO perturbative correction to r]b — ?■ J/'ipJ/tp J36|. 

V. PHENOMENOLOGY 

With our NLO calculations completed, following the formulas ([9]) and f[T3|) . we are ready 
to carry out a detailed analysis for the processes e'*"e^ — )■ J/ip^ijj') + Xcj and confront the B 
factory measurements. In our numerical analysis, we set ^/s = 10.58 GeV, m^ = 1.5 GeV. 
The electromagnetic fine structure constant is chosen as a{y/s) = 1/130.9 J39|. The running 
QCD strong coupling constant is evaluated by using the two-loop formula with A^ = 0.338 



GeV [111, |12[. The nonpertubative input parameters, i.e., the wave function at the origin 
for J/ip {ip'), as well as the ffist derivative of the radial wave function at the origin for Xcj, 
bear a fair amount of uncertainties. Their values have been compiled in Ref. [40|, which 
are estimated from several different potential models. We choose to use those given by the 
Buchmiiller-Tye potential model |4l|: |i?j/^(0)p = 0.81 GeV^ \R^'{0)\'^ = 0.529 GeV^ and 
\R' (0)P = 0.075 GeV^ 

' XcJ ^ ^ ' 

Another important source of uncertainty for the NLO predictions comes from the scale 
setting for the strong coupling constant. As is well known, this scale ambiguity is a typical 
nuisance of NRQCD factorization approach, reflecting the fact that two disparate hard scales, 
^/s and rric, are entangled together in NRQCD short-distance coefficients. In contrast, the 
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light-cone approach, when armed with the idea of refactorization, can efficiently disentangle 
these two scales, so one naturally expects that the scale ambiguity will be greatly reduced. 
While it may sound natural to choose the renormalization scale /i to be around the highest 
scale i/i, however, as is clearly illustrated in 18|, it is rather unappealing to set the scales 
entering all as in NLO short-distance coefficient to be around ^/s unanimously. The reason 
is that a part of NLO correction comes from the loop region with lower virtuality, and as such 
can be identified with the NLO perturbative correction to the charmonium decay constant, 
so the corresponding a^ should definitely be affiliated with a scale around 2mc rather than 

Since it is not possible to completely solve the scale ambiguity problem within the confines 
of the NRQCD factorization approach, we proceed to estimate the cross section by assign- 
ing all the occurring as with a common scale, /x, and choosing /z = i/i/2 and // = 2mc, 
respectively. It is hoped that the less biased results interpolate between these two sets of 
predictions. 



• JM+ztcO 



e*e ->J/i/f+^c 



• il<l'+X^i 




FIG. 5: The /i-dependence of LO and NLO cross sections for e+e^ — )■ J /il) + Xcj (J = 0, 1, 2) at 
^/s = 10.58 GeV. The uncertainty band is due to sliding rric between 1.4 GeV and 1.6 GeV, where 
the central curves correspond to rric = 1.5 GeV. 
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FIG. 6: The LO and NLO cross sections for e~^e~ — )• J/ip + Xcj (J = 0, 1,2) as a function of ^/s. 
The uncertainty band is obtained by varying ^ from ItUc to y^, where the central curves represent 
the default choice ^ = ^ . 



To develop a concrete feel about the scale dependence, in Fig. [5] we exphcitly show the 
/i-dependence of the LO and NLO total cross sections for the processes e~^e~ — )■ J/V' + XcJ 
(J = 0, 1, 2) at B factory energy. As can be clearly seen, including the NLO correction has 
notably reduced the scale dependence for J/ip + Xci,2, whereas of little impact for J/ip + Xco- 
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TABLE I: Polarized cross sections for each helicity channel and unpolarized (total) cross sections 
(in unit of fb) for e^e~ — )• J/ip + Xcj (J = 0,1,2). In the rightmost column, we also list the K 
factor for the unpolarized cross sections. We choose the foUowing input parameters: rric = 1.5 
GeV, fj. = a/s/2, accordingly, as(/x) = 0.211. 







O"(0,0) 


^(1,0) 


^(0,1) 


^(1,1) 


^(1,2) 


O"tot 


K 


J/i^ + Xco 


LO 
NLO 


1.11 
1.92 


1.86 
3.35 








4.83 
8.62 


1.79 


J/i' + Xci 


LO 
NLO 


— 


0.0012 
-0.0078 


0.37 
0.49 


0.033 
0.028 


— 


0.81 
1.01 


1.24 


J/i^ + Xc2 


LO 
NLO 


0.43 
0.44 


0.27 
0.33 


0.064 
0.081 


0.033 
0.039 


0.0023 
0.0026 


1.17 
1.33 


1.14 



TABLE II: The same as Table HI except we choose the different renormahzation scale: ^ = 2mc, 
accordingly Qs(^) = 0.259. 







O"{0,0) 


^(1,0) 


^(0,1) 


^(1,1) 


^(1,2) 


<7tot 


K 


J/i^ + Xco 


LO 
NLO 


1.67 
2.50 


2.80 
4.46 









7.26 
11.43 


1.57 


J/i' + Xc\ 


LO 
NLO 




0.0017 
-0.016 


0.56 
0.55 


0.050 
0.020 




1.22 
1.11 


0.91 


J/i) + Xc2 


LO 
NLO 


0.65 
0.40 


0.40 
0.35 


0.097 
0.089 


0.050 
0.041 


0.0035 
0.0026 


1.76 
1.38 


0.78 



Also in Fig. |5] we have examined the mc-dependence of the cross section, which is reflected 
in the error band by varying nic from 1.4 GeV and 1.6 GeV. 

In Fig. El we also plot the LO and NLO total cross sections for e^e~ — )■ J /ip + Xcj 
(J = 0, 1, 2) as a function of y/s. The error band is obtained by sliding /x from 2mc to ^/s, 
where the central curves represent the default choice jj, = ^. For the LO predictions, all 
the J /ip + XcJ channels have comparable widths of the error bands. But the NLO results 
for J /ip + Xc2 exhibit a significantly narrower band compared to those for J /ip + Xco,i- 

Table [T] and Table [TTl tabulate our predictions to both of the LO and NLO cross sections 
for e+e~ — )■ J /ip + Xcj {J = 0, 1, 2), with these two sets of /x chosen respectively. In addition 
to the unpolarized (total) cross sections, we also include the polarized cross sections from 
each helicity channel. 

Let us first discuss the unpolarized cross sections. For e"'"e~ — )■ J /ip + XcO; the NLO 
correction has significantly enhanced the LO prediction, with a K factor of 1.79 and 1.57 
respectively, corresponding to two different choices of ^ ^. In contrast, the NLO corrections 



^ We find disagreement with the previous NLO correction calculation for e+e -^ J ji]) + XcO 116|- When 
taking the same input parameters as theirs, we obtain the K factor of 1.57, while theirs is 2.8. 
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TABLE III: Comparison between our predicted unpolarized cross sections with the measurements 
at B factories for e+e^ — )■ J/ip{Tp') + Xcj (J = 0, 1, 2). The cross sections are in units of fb. We fix 
rric = 1.5 GeV, the error is estimated by varying /u from 2mc to y/s, where the central value refers 
to /x = -y/s/2. 





Belle 


BaBar 


LO prediction 


NLO prediction 




axe>2(o)[42] 


a X %2[6] 






(^{J/'ip + Xco) 


6.4 ±1.7 ±1.0 


10.3 ±2.5+1-8 


A SQ+^-43 


8.62+^-^0 


<J/i^ + Xci) 


- 


- 


n O1+0.41 


1 m+o-io 


(y{J/i^ + Xc2) 


- 


- 


1 17+0.59 
■•-•^'-0.38 


-, qq + 0.04 

i.OO_Q20 


<y{.JI^ + Xci) + (y{Jli> + Xci) 


< 5.3 at 90% CL 


- 


1 QS+lOO 


9 qr:+0.14 
^••-'^-0.38 


(y{^' + Xco) 


12.5 ±3.8 ±3.1 


- 


2 79+^-41 


4.98li:i 


f^(^' + Xcl) 


- 


- 


47+0-24 


0.58lH^ 


C^(^' + Xci) 


- 


- 


n co+0.34 
0.68_o22 


0.77in^ 


(tO^' + Xcl) + (Jii^' + Xc2) 


< 8.6 at 90% CL 


- 


1 14+0.58 


1.36+°-°« 



to e+e" — )■ J/ip + Xcl, 2 have a milder impact, even with the sign uncertain. Concretely 
speaking, the first setting of fi tends to increase the LO results modestly, while the second 
setting tends to reduce the LO results to some extent. This behavior may be clearly seen in 
Fig.O There the LO and NLO prediction bands cross with each other for e+e" — ?■ J/'ip+Xci,2, 
which implies that the K factor could be above or below 1, depending on the chosen scale. 
A similar understanding can also be achieved in Fig. [61 At ^/s = 10.58 GeV, one sees 
from Fig. [6] the uncertainty band of the NLO predictions for e+e^ — )■ J/ijj + Xci,2 has been 
completely submerged inside the band of the LO predictions. Therefore, depending on 
whether choosing the lower or upper bound for LO cross section, the K factor would be 
greater or less than 1. 

In Table UTTl we compare our predicted cross sections with the measurements at B factories 
for e+e~ — )■ J/ipi'ip') + Xcj {J = 0, 1,2). One sees that the NLO perturbative correction to 
J I^J + XcQ really brings the NRQCD prediction closer to the data, although there exists some 
slight tension between BaBar and Belle measurements. By far, the e+e" — > J /ip{ip') + 
Xci,2 processes have not yet been observed in any experiments. Our predictions for the 
J/ip{ip') + Xcl, 2 cross sections are compatible with the upper bounds placed by the Belle 
experiment. However, even if the large positive NLO correction is taken into account, our 
predicted ip' + Xco cross section is still significantly below the central value of the Belle 
measurement. To clarify this puzzling situation, it seems necessary, and, urgent, for BaBar 
to perform an independent measurement for this process to see whether it confirms or 
disconfirms the Belle results. 

It seems foreseeable that, in future Super B factory, experimentalists may be able to 
measure some of the polarized cross sections for the e+e~ — )■ J/ip + XcJ processes. Therefore, 
it is informative to examine the polarized cross sections in Table |T] and Table HTl 

For J/ip + Xco production cross section, both the (0, 0) and (1, 0) hehcity channels have 
comparable magnitude, either at LO or at NLO, and the latter is even somewhat greater. 
This is quite counterintuitive, diametrically contradicting what is expected from the helicity 
selection rule. This might be viewed as a hint that the B factory energy may lie still far 
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from the asymptotic scaling regime. 

For e"'"e~ — )■ J /ip + Xci, the hierarchy among the cross sections of three different hehcity 
channels is also somewhat abnormal. The contribution from the (0, ±1) channel is far more 
significant than the other two, and to a good approximation, one only needs to retain 
this channel. This pattern still holds true after incorporating the NLO correction for both 
settings of /i. It will be interesting for future Super B experiments to verify that the angular 
distribution of J /ip (or Xci) is predominantly of form 1 + cos^ 9. In passing, we also note 
that the tiny LO contribution from the (±1,0) state may be attributable to the accidental 
suppression factor with respect to the helicity selection rule. After including the NLO 
correction, cr(i,o) even becomes negative. This can be easily understood by inspecting (I15al) . 
since the NLO helicity amplitude is 1/r enhanced relative to the LO one, and with opposite 
sign. However, in practice this polarized cross section is too small to be measured. 

For e^e~ — )■ J /ip + Xc2i the hierarchy among the cross sections of five different helicity 
channels roughly obeys the helicity selection rule, except the contribution from the (0, ±1) 
channel is much smaller than the (±1,0) channel. From Table [I] and Table [TTl one sees that 
the bulk of the cross sections comes from only two hehcity states, i.e. (0,0) and (±1,0). 
For IX = y/s/2, the NLO correction has small impact on both helicity states; for /x = 2mc, 
the NLO corrections push down both polarized cross sections to some extent. It would be 
interesting for the future high-statistics experiment to observe this production process, and 
test whether the produced Xc2 are predominantly longitudinally-polarized. 



VI. SUMMARY 

In this work we have computed the complete NLO perturbative corrections to e"^e" — t- 
J/i^ + Xcj {J = 0, 1, 2) within the NRQCD factorization framework. The NLO NRQCD 
short- distance coefficients are inferred by directly extracting the contribution from the hard 
loop-momentum region. We have calculated the NLO corrections to each of the 10 inde- 
pendent helicity amplitudes. We have made a detailed analysis for both the polarized and 
unpolarized cross sections and compared with the measurements at B factory. 

We find a significant positive NLO perturbative correction to e"'"e^ — )■ J/ip ± XcOy which 
helps to bring the predicted cross section in agreement with the B factory measurements. In 
contrast, our NLO predictions to ip' ± Xco cross section is still significantly below the central 
value of the Belle measurement. It is still too early to draw any solid conclusion about 
whether NRQCD factorization fails for this channel or not. We perhaps need to wait until 
BaBar collaboration carries out an independent measurement for this process. 

The impact of NLO corrections to e+e~ — >■ J/ip + Xci,2 seems to be rather mild, even 
with their signs uncertain. Notice the predicted cross sections for these processes are about 
5 or 6 times smaller than that for e+e^ — t- J/tp ± Xco- Hopefully the future Super B factory, 
with much higher luminosity, will eventually observe these two channels. 

Our studies of polarized cross sections reveal that the bulk of the total cross section comes 
from the (0, ±1) helicity state for e~^e~ — )■ J/ip + Xci, and from (0,0) and (±1,0) helicity 
states for e+e~ — )• J/ip + Xc2- It will be interesting for the future Super B experiments to 
test these polarization patterns. 

On the theoretical side, we have worked out the explicit asymptotic expressions of all 
the ten NLO hehcity amplitudes for the e+e~ — )■ J/ip + Xcj {J = 0,1,2) processes. We 
confirm that helicity selection rule is modified logarithmically at NLO in ag- The pattern 
we recognize in these NLO asymptotic expressions lends further support for the speculation 
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made in Ref. [18[: The hard exclusive processes involving double- quarkonium at leading twist 
can only host the single collinear logarithm In s/m^ at one- loop, while those beginning with 
the higher twist contributions are always plagued with double logarithms of form In^ sjw?^. 
It is of some theoretical interest to reproduce the asymptotic expressions for the helicity- 
conserving channels such as e+e^ — )■ J/'?/'(0) + Xco,2(0) in the light-cone approach. This 
should be definitely feasible, which is guaranteed by the leading-twist factorization theorem. 
Nevertheless, it is much more challenging for the light-cone approach to reproduce, and 
resum, those pro cess- dependent double logarithms associated with the helicity-suppressed 
channels. 

ISiote added. After the calculation was finished and while we were preparing the draft, a 
related work appeared in arXiv recently [43|, which also investigated the 0{as) correction 
to the processes e~^e~ — )■ J/ip + Xcj [J = 0, 1, 2) within NRQCD factorization approach. 
Once taking the same input parameters as theirs, and summing up the contributions from 
all the helicity states, our results agree with Ref. |43| on the numerical sizes of the NLO 
corrections to the unpolarized cross sections for each J = 0, 1, 2. 
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